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We consider a very large ultra-relativistic nucleus. Assuming a simple model of the nucleus and 
weak coupling we find a classical solution for the gluon field of the nucleus and construct the two- 
dimensional color charge density for McLerran-Venugopalan model out of it. We prove that the 
density of states distribution, as a function of color charge density, is Gaussian, confirming the 
\Q ' assumption made by McLerran and Venugopalan. 
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(]J 1 I. INTRODUCTION. 

on : 

, Consider a very large nucleus, probably larger than can be physically realized. The nucleons are distributed 
homogeneously inside the nucleus. Recently L. McLerran and R. Venugopalan proposed a program of computing the 
s ! | gluon distribution fuctions for such a nucleus at small x [1,2]. 

One of the interesting problems in the McLerran-Venugopalan model for the small- a; part of the gluon distribution 
of such a large nucleus [1,2] is finding the classical solution for the gluon field treating the valence quarks of the 
nucleons in the nucleus as recoiless sources, which are delta functions along the light-cone when the nucleus is moving 
\f} ■ near the velocity of light. A convenient way to deal with the problem is by working in the light-cone gauge. The 
source is characterized by a two-dimensional color charge density p{x), where x is a vector in the transverse direction. 
The proposed model assumes that in order to find the average value of any observable having longitudinal coherence 
length long compared to the nucleus, one calculates this observable for a given p(x), and then averages it over all p 
^ ■ with the measure 

J [dp] cxp (~ I d 2 xp 2 (x)^j , (1) 

. . , where p is the average charge density squared. 

We consider a large nucleus consisting of "nucleons" , which for simplicity of description are chosen to be just 
quark-antiquark pairs (see Fig. 1). Valence quark and antiquark are treated as point particles free to move inside of 
$_j ■ the nucleon, but unable to get out. 
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FIG. 1. Nucleus with "nucleons" being quark-antiquark pairs. 



We are interested in the gluon field of the ultrarelativistic nucleus viewed in the laboratory frame. We assume 
that the field in each individual nucleon is not large. This allows us to approximate the covariant gauge potential of 
each quark by a single gluon exchange. In a covariant gauge the classical field of a single ultrarelativistic particle is 
proportional to a delta- function in the X- direction [3]. Since in our model of ultrarelativistic nucleus different quarks 
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have different X- coordinates in the laboratory frame, the fields of individual quarks don't overlap, which allows us 
to superimpose them and justifies our single gluon exchange approximation. Then the total field of the nucleus in the 
covariant gauge is the sum of the quark fields. We make a gauge transformation which changes the total potential to 
the light-cone gauge. So, we get a solution to the classical equations of motion in light-cone gauge, where the field, 
A lll is directly related to the gluon distribution in the small- a; region [1]. 

Following McLerran and Venugopalan the nucleus is considered to be very large, thus although the field of each 
individual nucleon is weak, the total field is strong at low momentum in the light-cone gauge due to the overlap of 
the fields of a huge number of nucleons. Still we can neglect the contributions of several nucleons, without changing 
the answer, i.e., we work in the leading power of the number of nucleons. 

We may treat the source as classical only when we are at sufficiently small momenta that the individual quarks can 
not be resolved. It was shown in [1] that this requires that k 2 -C p 2 , where k is the typical momentum scale. The weak 
coupling approximation is valid when fc 2 3> ^qcd- Then the momentum range we consider is ^qqd <^ k 2 p 2 . 

Now the task is to construct the two-dimensional charge density, giving the correct classical solution. This is done 
by just substituting the classical solution in the equations of motion. The density we find this way happens to satisfy 
the Gaussian distribution. We show this by calculating the correlation functions of the densities at different transverse 
points and proving that they are exactly what one would expect for the Gaussian distribution (^). That is we justify 
the method for averaging the observables proposed by McLerran and Venugopalan. 

In sect. 2 we calculate the solution of the classical equations of motion - the non-Abelian Weizsacker- Williams field. 

In sect. 3 we construct the two-dimensional charge density. 

In sect. 4 we show that the charge density has a Gaussian distribution by calculating the correlation functions. 
In sect. 5 we confirm the techniques proposed in [1]. 



II. APPROXIMATE SOLUTION. 



We start with some random distribution of nucleons in the nucleus and quarks and antiquarks in the nucleons. 
The nucleons , in the rest frame of the nucleus, are assumed to be spheres with equal radius, and the quarks and 
antiquarks are distributed randomly inside each sphere, with equal probability to be at any place inside the sphere, 
but with zero probability to get outside. The density in the rest frame is given by 

8 

P (f) = ^r/(f) (2) 

o=l 

with 

N 

P a (x) =gJ2(T?)[6(x--x i ) 5{x - %)], (3) 

z=l 

where af{ is the coordinate of a quark in the ith nucleon (there are N nucleons in the nucleus), x\ is the coordinate 
of the antiquark, T a are generators of SU(3) in color space, (T?) are similar generators in the color space of each 
nucleon. The reason we separate them is because p a comes from the current = gq a ^^(T a ) a (jqp, so the expression 
for p a should include a T a acting in each individual nucleon's color space. 
In the laboratory frame for the ultra-relativistic nucleus the density is 

8 N 

p(x,x^) = ~%J2J2 T a ( T i)[ s ( x - - x -i) s ( x - x i) - 6 i x - - x'_i)8{x - o£)]. (4) 
V z a= i i= i 

Assuming that the coupling is weak and using the expression for the potential of a single particle in a covariant gauge 
[3, Appendix A] we approximate the field of the nucleus by superposition as 

8 N 

-8(x--x'_ i )hi(\x-£\\)] ,A' = Q,A'_ = 0, (5) 

where A is some infrared cutoff. The prime at the field or the field strength denotes the covariant gauge. The field 
strength in the covariant gauge is then 
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8 N 

K± = i E E ra ^) ( - ^Sfe - - ^iffffe ) ■ (6) 

a— 1 i— 1 



From here on the subscript _L will mean that the object is a vector in the transverse space over this index. 

We now do a gauge transformation to transorm this field into the light-cone gauge. The potential in a new gauge is 

= SA'^S- 1 - ^(d^S)S- 1 . (7) 

Requiring the new gauge to be the light-cone gauge, A + = 0, we obtain 

S(x, x_) = P exp ^— ig J dx'_A' + (x, x'_)^j . (8) 

Then the field in the light cone gauge is 

/X- rX- 
dx'_F +J _{x,x'_) = / dx'_S{x,x'_)F' +x {x,x'_)S~ l {x,x'_), (9) 
-oo J — oo 

Only transverse components are non-zero. Substituting F[ j_ (x, x'_ ) from (0) we get 



8 N 

A{x,x.) = |^EE^ Q ) ( S(x, x^TS-^x, x.i) 6{x - ~ 

a — 1 i—1 



(10) 



This is our estimate of the solution of the classical equations of motion for a given configuration of quarks inside the 
nucleons and nucleons inside the nucleus. Formula gives us the non-Abelian Weizsacker- Williams field generated 
by the valence quarks. 

III. TWO-DIMENSIONAL COLOR CHARGE DENSITY. 

The equation of motion is 

D^F^ = J v . (11) 

In the McLerran-Venugopalan model [1] the classical current J M only has components in the + direction and is 
proportional to a delta function of x_ : 

J*{x) = 5» + 5{x„)p(x). (12) 

This can be treated as a definition of the two-dimensional color density. Our goal now is to construct p(x). Integrating 
both sides of ([fl]) over x_ and using (|l^) gives 

/ + OC 
dx-DiF +i (x,X-) 
-oo 

f + OO 

dx-{diF + i(x,X-) - ig[Ai(x,x-),F+i(x,x-)]}, (13) 

where i = 1,2 (transverse direction). Using A(x, a;_) from ( [To|) , we can calculate F + ±(x, X-). Substituting both in 
( |l3| ) we end up with the following expression for the density: 

8 N 



P(x) = .9 E E(^)[ 5 fe' x-JT'S- 1 ^, x-i)6(x - x t ) 

a—1 i—1 

-S(^,x'_ i )T a S- l ^,x , _ i )S(x-^)]. (14) 
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The details of calculations are presented in Appendix B. We can see now that our expression for two-dimensional 
density is just a rotation of the three-dimensional density (in the laboratory frame) we started with: 

/+oo _ 
dx-S(x, X-)v2p(x, a;_)5 _1 (a;, xJ). (15) 
-oo 

In the expression for the light-cone potential we had an infrared cutoff A, so it may seem that this cutoff will appear 
in the S(x, and, consequently in p(x). However this is not the case, because, although the light-cone potential 
is cutoff dependent, S(x,X-) is not. To see this let's do an explicit calculation of S(x,x_): substituting (||) into (||) 
and using the definition of the path-ordered exponential we obtain 

N 

S(x,x-) = Y[ (0(x'_i - x-i)e^e^ + 9{x_ t - .xLje^e^) , (16) 

i=l 

with 

E » = ^E Ta (^ )M|z-^|A)%- -x-i), 

Zn * — ; 



(17) 



But the matrices we exponentiate (E, and E£) commute (for the same i), so 



iV 



S(x, x-) — JJexp 



2tt 



a=l 



k-Xi\ 



(18) 



Here we neglected the contribution of the "last" nucleon, i.e. the nucleon (or several nucleons) whose quarks or 
antiquarks may overlap the point x_ at which we calculate S(x, X-). These nucleons may potentially cause us some 
trouble, but the philosophy of the large nucleus approximation implies that the fields of individual nucleons are small, 
and we construct a strong field out of a large number of nucleons. The contribution of each individual nucleon is 
negligible, it's their sum which matters. That means we can neglect these "last" nucleons. 

Another way to say this is that we want to do a calculation keeping the leading powers of N only. Then dropping 
few nucleons will not change our result. 

From ( |l8| ) we see that S(x,X-) is cutoff independent, and so is the density p(x). 



IV. CALCULATION OF CORRELATION FUNCTIONS. 



Now that we found the charge density, let's show that it's distribution is Gaussian by calculating the density 
correlation functions. First we note that the average density is zero, as expected: < p a (x) >= 0, where < . . . > 
denote the averaging over all possible positions of quarks and antiquarks in the nucleons, and nucleons in the nucleus, 
as well as averaging over all possible colors (keeping each nucleon color neutral). 

For two densities correlation function we have : 

< p'(^(»j >= n / T^i0?^^mm. (M) 

where R is the radius of the nucleus, a is the radius of the nucleons, is the position of the center of the fcth nucleon 
in the nucleus (in the rest frame), Xk and x' k are positions of the quarks in the nucleons; (aa\ . . . \/3f3) implies an 
average over all color-neutral states of the nucleons. 
Using @ we obtain (since p a (x) = 2Tr[T>(x)]) 
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x {2Tr[T a S(x i ,X- i )T c S- 1 (x i ,X- i )]S(x-x i ) - 2Tr[T a S(x! i ,x'_ i )T c S- 1 (x! i ,x'_ i )]6(x-x! i )} 



x {2Tr[T b S(x j ,X- j )T d S- 1 (x j ,x^ j )}S(y - % ) - 2Tr[T b S(x! j ,x'_ j )T d S- 1 (x! j ,x'_ j )}6(y - x^)} \(30). (20) 

In S(x^x-i) the "last" nucleon is the ith nucleon. Applying the same arguments we had before we can drop this 
"last" nucleon. Then there will be no (T?) matrices in and S(gc^,x'_ i ). It is convenient to label nucleons 

according to their coordinates along the X- axis. The greater the coordinate, the greater the number of the nucleon. 
Without any loss of generality we can assume that i > j. Then if i > j the (Tf) matrix in front will be the only 
matrix in the ith nucleon color space for the term corresponding to fixed i and j in our expression. So, when we do 
the color averaging it will give zero [Tr(Tf) = 0]. That means that i — j. Then color averaging in the ith nucleon 
space gives ^-TT[(T[)(Tf)} = ^-S cd . 

Each density is a difference of the quark and antiquark parts. In ( po| ) one can easily see that the product of quark 
(and antiquark) components gives a delta-function of x and y , while the product of quark and antiquark components 
gives some smooth function. In most of the physical applications we'll be looking at the scales much smaller than the 
nucleon's radius a [4]. At these scales we can neglect this cross terms with respect to the delta- function terms: 

X ' d 3 r t ^3 r {x- Li y ^ f d 3 x k d 3 x' k 



/ — 1 i — 1 K, — X 



[(4/3)7ra 3 ] 2 

x(aa|^4Tr[r a 5(x l ,x_0r c 5- 1 ( S ,^0]Tr[r b 5( S ,x_0r c 5- 1 ( S ,a;_0]|/3^). (21) 

c=l 

After dropping the "last" nucleon S(x, X-i) = SiXjX 1 ^). 

For two traceless 3x3 matrices M and N the following formula is true: 
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Tr [MT a ] Tr [N T a ] = -Tt[MN}. (22) 



a=l 



Using this, and making some approximation when integrating over fj (a <C R), we get 



<P a ^ P Hy)>=S a %^Six-y) ]/ l-^. (23) 

Here , for the first time we made an assumption about the geometry of the nucleus and nucleons in it - when doing 
the average over the positions of quarks and nucleons we assumed that the nucleus and nucleons are spherical in the 
rest frame. For a cylindrical nucleus (in ^-direction) one would get < p a (x)p b (y) >= S ab jj-^z-5(x — y). 

By employing a similar techniques of dropping the "last" nucleon in the Six^jX-i) and keeping the leading powers 
of N only, we can show that the four-density correlation function is 

< P a (x)p b (y)p c (z)p d (w) >=< P a {x)p b {y) X p c (z)p d (w) > + < p a (x)p c (z) > 

x < p b {y)p d {w) > + < p a {x)p d {w) >< p b (y)p c {z) > (24) 

and prove a similar formula for a correlation function with any even number of densities, i.e., that it can be represented 
as 

< p ai (x})...p a2 -(x 2n ) >=<p ai (x 1 )p a2 (x 2 ) > ... < p 02 ™- 1 ^ 2 ™- 1 )// 12 ' 1 ^ 2 ™) > 

+ permutations, (25) 

where x , . . . ,x 2n are just 2n arbitrary points in the nucleus. Also we can show that a correlation function with an 
odd number of densities is zero to all orders in N (for details see Appendix C): 

< P ai (x 1 ) ■ ■ ■ P a2n+1 {x 2n+1 ) >= 0. (26) 
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V. CONCLUSIONS. 



By calculating the two densities correlation function (E3J) and proving (25) and (Eq) we showed that for any n points 



<p ai fe 1 )...p a " (£")>= 7 \ MX (27) 

/[dp] exp (^/^^g. 



with 



M l - J 2 N c ttR 2 \ R 2 2 iV c ™ 2 V R 2 [ ' 

in our model. Note that our p 2 goes as TV 1 / 3 , as expected. The average on the left hand side is understood as 
< . . . >= nit=i / (a%)kR* [(4/3)^0^^ ( a ® \ ■ ■ ■ IftP)- F° r a spherical nucleus p 2 is a function of x. 

The most general observable in our system can be represented as some functional of J K{x' ,x)p{x)d 2 x, namely the 
value of this observable for a given p(x) is 

^■'(/"'H' (29) 

where F(f) is some functional and K(x',x) is some kernel independent of p(x). For instance, one can consider 
Wilson loop: W =< Tr Pexp(— ig J dx • A) >. It's a functional of A(x). The field A(x) can be represented as 
J K(x' , x)p(x)d 2 x. Then, using the definition of a path-ordered integral, we can expand the Wilson loop in the 
powers of A(x): 

W =<Ty]\[1 - igdx k ■ A{x k )] >=<Tr^ Cn A(xi)...A(x n ) >, (30) 

k n 

with some coefficients c„. Writing A(x) as A(x) = J2a=i T a J K a (x,x')p(x')d 2 x' , we achieve 

8 

W = Y^c n T4T ai ■■■T a "} d 2 x[...d\l l K a i(x 1 ,^ 1 )...K a "(x n ,x! n )x 

x <p(x! 1 )...p(x! n )>. (31) 



Now we can use (g7|) , getting 



/[dp] exp (- / d 2 x^gj) Tr P exp (-ig J dx • A) 
^" J[dp]exp(-/^|g) ' ^ 



A similar treatment can be applied to any observable ( which can be represented as (|29|)) to prove that 

f[dp] exp 

This confirms the assumption made by McLerran and Venugopalan in their model [1] Q 



*I have recently learned (private communication from L. McLerran to A. Mueller) that similar conclusions have been reached 
by Professor L. McLerran and collaborators. 
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APPENDIX A: 

From [3] the light-cone potential of a charge e moving with a velocity v is 

A + = ^ 1 ,A = A_ = 0. (Al) 

47r hx 2 _ + (1 - v 2 )x 2 



If we do a Fourier transform 



A„(k) = y-^rz I dx + dx-d 2 xe lk + x - +lk - x +- l ^A^x) (A2) 



(2tt)* 

and take the limit v —> 1, we get for the ultrarelativistic particle: 



A + (k) = & -^>. (A3) 

Z7TK 



If we go back to the coordinate space, by doing an inverse Fourier transform of (A3), we end up with 



A + (x) = -^-S(x^)\n(\x\\), (A4) 
which is different form the v —> 1 limit of the original expression by a gauge transformation. 

APPENDIX B: 

Defining 

8 

fi(x,x-) = ±J2(7?)S(x,x-)T a S- 1 (x,x-) (Bl) 

Z7T * — ' 

a=l 

we can rewrite ( fI0| ) as 

N / - - 1 \ 

^ V m -x.il 2 k-x'^ 2 j 

Then the field strength is 

d N ( x-x 

F +± (x,X-) = -^A{x,x-) = X ( fi(x,x-i) ~_ ~\ 2 S ( X - ~ x -i) 



,\A X - x '-i) ■ (B3) 



Substituting QBJ) and (jBJ) into (|3|) we get 
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p{x) = 'Y^[fi(x,x^ i )2'K8{x-x j ) - fi(x,x'_i)2TrS(x - 
N / - ■ 



i=i 

N 



~ { 9 X [[fi(£> x -i)ifj(M-> x -j)h 



\x-x^\x-xA 2 " [X - j X ~ l > 
- [/i (a ar-i), /j ix, x' )} ,f ~| 2I ~ "^ ^-j ~ 
■[/jfe ar_j)i /jfe^-i)] ^ _ p _ ^2^( x -i _ x -i) 



+[f i (x, x '_ i ),f j ^ x '- j )} ] :_: / ) 2 ] ^_% o(x'_ j -x'_ i ) ) . (b-i, 



A straightforward calculation yields 

3 8 8 JV 



X — Xa 



and 



where 



v ' 0=1 6=1 j=l 1 

)T b W(x,X-i,X-j )T a S f (a;, — S(x,x—i)T a W(x,x-j,x-i) 
x r^fez-i)] - .-~-.\M x-i - x'_j)lS(x, x'_ j )T b W{x,X- i ,x'_ j )T a 

\X — Xj\ J J J 

x S^&x-i) - S(x,X- i )T a W(x,x'_ j ,x- i )T b S- 1 (x,x'_ j )]} (B5) 



2 8 8 

Ifi&x-ilfs&x-j)] = 1 ^^^(Tt)(Tj)[S^,x- i )T a W(x i x-j,x- i ) 

\ ' o=l 6=1 

x T 6 S -1 (x,a;_j) - S(z, sg_.,-)T 6 W(x, z-i, x-j)T a S' l {x, x-i)], (B6) 



W(ac,a;_i,a;_j) = Pexp — ig / efe_A+(x, x_) . (B7) 



Plugging this back into (B4) we end up with 

AT 

P(x) = X^fe ^-02^fe - ^ - /ife x'_,)2 7 r5(a I - ^)], (B8) 

i=l 

which is equivalent to (|l4j). 

APPENDIX C: 

Let's prove that the three-densities correlation function is zero. Then the proof for an arbitrary correlation function 
of an odd number of densities will become obvious. Similar to (pet) we write 
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N f j3 j3 j3 / 8 w 

< /fey (jKb) >- f n / T3OTTWW ( " S| S S roWx* > 

i=l' / w y LW / J o',6',c'=l*,j,fc=l 

x (2Tr[T a S{x tl x-JT*' S^fa, X-i)]5(x - S ) - 2Tr[T a S&, x'^' S~\^ *'_*)]% - a£)) 

x (2Tr[T b 5(x J ,^)T b '5- 1 (x J , a ;_ J )]% -^.) - 2Tr[T fc 5(^, x'_ )T b ' S' 1 ^, x'_ )]5{y - a£)) 

x (2Tr[T c ,S( ; r fe , a ;_ fe )r c ^- 1 (x fc , a ;_ fc )] ( 5(z -a*) - 2Tr[T c S(^, a^T^S" 1 ^ a/_ fc )]«U - a&)) » (CI) 
After dropping the "last" nucleon and averaging over the colors of this nucleon, we get 

i ri d 3 xid 3 x[ 



< 



d 3 ri d 3 xid 3 3 



fJ-J (4/3)7ri? 3 [(4/3)7ra 3 ] 2 
x (aa\ ^^[{Tf){T?){T?)XTr[T a S x T a ' ' S- 1 ]Tr[T b S y T b ' S-^Tr^S^ ' S' 1 }^) 



8 N 



a' ,b' ,c' — 1 i— 1 

x [S(x xj S(x a£)] [% - s.) - 6(y ^)]\6{z - x*) 6(z (C2) 

where S x = S(x, X-i) = S^x'^) (after dropping the "last" nucleon) and independent of x±. The product of three 
brackets with delta-functions integrated over x^ and x^ obviously gives zero. So, < p a (x) p b (y) p c (z) >— , as 
advertised. 

Similar techniques can be applied to an arbitrary odd number of densities to show that their correlation function 
is zero. 
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